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Abstract 


The performance of a direct alkaline fuel cell (AFC) is studied separately using methanol, ethanol and sodium borohydride as fuel. Potassium 
hydroxide solution was used as an electrolyte. Pt-black and manganese dioxide catalyst were used to prepare the anode and cathode electrodes. 
Ni mesh was used as current collector. The direct alkaline fuel cell was constructed with the prepared anode and cathode electrodes and Ni mesh. 
The current density—cell voltage characteristics of the fuel cell were determined by varying load and at different experimental conditions, e.g., 
electrolyte concentration, fuel concentration and temperature. The fuel cell performance increases initially with the increase in electrolyte (KOH) 
concentration and then decreases with further increase of the same. The cell performance increases initially and then no appreciable improvement 
noticed with the increase in fuel concentration. The performance of the fuel cell increases with increase in temperature in general with the exception 
to NaBH; alkaline fuel cell. A mathematical model for the direct alkaline fuel cell is developed based on reaction mechanism available in the 
literature to predict the cell voltage at a given current density. The model takes into account activation, ohmic, concentration overpotentials and 
other losses. The model prediction is in fair agreement with the experimental data on current-voltage characteristics and captures the influence of 


different experimental conditions on current-voltage characteristics. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Alkaline fuel cell (AFC) was the first fuel cell technology 
used in many practical applications like Apollo space missions 
and automobiles [1,2]. The development of AFC technology 
had reached its peak in the beginning of 1980s [3] but its further 
development was stopped due to many technical, commercial 
and safety issues. In last few years, the interest on AFC increased 
due to more favorable oxygen reduction [4,5] and fuel oxida- 
tion reactions [6] in alkaline condition. Apart from these, the 
cost, simplicity, efficiency and the possibility of use of non- 
noble metal catalyst [3,7] compared to other low temperature 
fuel cell technology have given an impetus to the AFC research. 
The detailed comparisons of AFC and polymer electrolyte mem- 
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brane fuel cell (PEMFC) technologies based on hydrogen fuel 
is given by McLean et al. [8]. Though, the hydrogen as a fuel 
has many advantages but the certain obstacles have led to look 
into possibility of direct use of hydrogen-rich liquid and solid 
fuels in fuel cells. Some of the obstacles associated with the 
use of hydrogen as fuel are the generation of pure hydrogen in 
large quantity at a lower cost comparable to gasoline cost, on 
board storage of hydrogen (low density), low power output per 
unit weight of the fuel cell and fuel processor and safety issues 
[2]. The investigators are working on direct feeding of hydrogen 
rich fuels, e.g., alcohols, ethers, N2H4, NH3, and NaBH; into 
the alkaline fuel cells [9—14]. These fuel cells are operated at 
a low temperature (20-60 °C) and produce low power density. 
Thus, direct alkaline fuel cell is suitable for use as power source 
to portable electronic equipments. 

Among the fuels, methanol is an attractive liquid fuel because 
itis relatively cheap, readily available, solubility in aqueous elec- 
trolytes, easily stored and handled. Energy density of methanol is 
6kWhkg~! (5 kWh L7 !). Ideally, the electrochemical oxidation 
of methanol produces six electrons per molecule of methanol. 
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Nomenclature 

A area of the electrode m?) 

C concentration of oxygen at cathode catalyst sur- 
face (kmol m~?) 

CE concentration of ethanol (mol L7!) 

CM concentration of methanol (mol L7!) 

Co concentration of oxygen in air (kmol m~*) 

COH concentration of KOH (g cc!) 

Cou concentration of KOH (mol L7!) 

Cep concentration of sodium borohydride (mol LT!) 

Deft effective diffusivity of oxygen from electrode 
composite (m? s7!) 

DA air diffusivity of oxygen from air (m? s7!) 

E open circuit voltage (V) 

Ecel cell voltage at any load (V) 

F Faraday constant (C mol`!) 

J current density (A m`?) 

Jo exchange current density (A m~?) 

kj rate constant of forward reaction in Eqs. 
(16)-(18),  (23)-(25) and (27)-(29), i=1-9 
(Lm~?s7!) 

Kg constant for ethanol oxidation in Eq. (26) 
(CL? mol-? m~? s7!) 

KM constant for methanol oxidation in Eq. (21) 
(CL? mol-? m~? s7!) 

Ksg constant for sodium borohydride oxidation in Eq. 
(30) (CL? mol~? m~? s7!) 

Ki adsorption-desorption equilibrium constant in Eq. 
(22) (i=1, 2), 26 (i=4, 5) and 30 (i=7, 8) 
(L mol!) 

k; rate constant of backward reaction in Eqs. (16) 
and (17), (23)-(24) and (27)-(28), i= 1-8 G43, 
6) (mol m~? s7!) 

k? standard rate constant in Eqs. (19), (26) and (30) 
(i=3, 6 and 9) (mol m~? s7!) 

L parameter in Eq. (A.1.1) (V) 

M constant, M = AbCo/n Fe?!" Do air 

n number of electrons 

N parameter in Eq. (A.1.1) (dimensionless) 

No oxygen molar flux (kmol m~? s7!) 

ri rate of the reaction in Eqs. (16)-(18), (23)-(25) 
and (27)-(29) (i= 1-9) (mol m=? s7!) 

R gas constant (J mol`! K7!) 

Ron, RA area specific resistance (Q cm?) 

T temperature (K) 

Tı temperature (°C) 


Greek letters 


a 


äm Ge 


transfer coefficient for the oxidation of different 
fuels (methanol, ethanol and sodium borohydride) 
transfer coefficient for oxygen reduction 
thickness of cathode (m) 

porosity of cathode (dimensionless) 
overpotential (V) 

activation overpotential (V) 


Neonc concentration overpotential (V) 


NFL undesired losses (V) 

Noh ohmic overpotential (V) 

Ocm,co fractional coverage of adsorbed CH3CO species 
(dimensionless) 

@yB0, fractional coverage of adsorbed HBO 2 species 
(dimensionless) 

Onco fractional coverage of adsorbed HCO species 
(dimensionless) 

Don fractional coverage of adsorbed OH species 
(dimensionless) 


However, methanol has disadvantages such as high toxicity, low 
boiling point (65°C) and not a primary fuel. Therefore other 
alcohols, particularly those derived from biomass resources are 
considered, as they are renewable in nature. Ethanol is one of 
the potential fuels as it is easily produced in large quantity by 
fermentation of biomass. The energy density of ethanol is about 
7.44kWhkg7! (5.9 kWh L7!) and it produces 12 electrons per 
molecule of electro-oxidation of ethanol. It is a non-toxic chem- 
ical and its boiling point is relatively higher compared to that 
of methanol. Ethanol and the higher alcohols have the draw- 
back that their C-C bond rupture is difficult in the presence 
of noble metal catalyst in the temperature range of 25—90°C 
[15]. Among solids fuels, sodium borohydride is gaining atten- 
tion because it contains 10.6-wt.% hydrogen and 7 kWh g7! of 
energy density, which is much more than the most of the hydro- 
gen storage alloys. The electro-oxidation reaction products of 
NaBH; are environmentally safe. Although hydrogen as a fuel is 
easily obtained from the hydrolysis reaction of borohydrides, the 
direct anodic oxidation of borohydride provides more negative 
potential than that of hydrogen gas. From the above discussion, 
it is clear that methanol, ethanol and sodium borohydride has 
excellent potential as fuel for direct alkaline fuel cell. 

In the present investigation a direct alkaline fuel cell is con- 
structed and its performance in the form of current—voltage 
characteristics is determined at different experimental condi- 
tions, e.g., fuel and electrolyte concentrations and temperature. 
The modeling of direct alkaline fuel cell based on methanol or 
ethanol or sodium borohydride as fuel is not reported in the open 
literature as far as the knowledge of the investigators. A math- 
ematical model on direct alkaline fuel cell is developed and the 
model predictions are compared with the experimental data on 
current-voltage characteristics at different operating conditions. 


2. Experimental 
2.1. Materials 


The catalyst used for the preparation of anode, Pt-black, was 
procured from Johnson Matthey Inc., UK. Manganese diox- 
ide used in the cathode was procured from E. Merck. The 
carbon paper (484C-1, Lydall, USA) was used as a substrate 
for catalyst particles. Nafion® (SE-5112) and Teflon® disper- 
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sions used as adhesives were obtained from DuPont, USA. Pure 
nickel screen (mesh 10) was used as current collector because 
it is highly resistive to alkaline condition and good electrical 
conductor. Analytical grade of methanol, ethanol and sodium 
borohydride were obtained from E. Merck. Potassium hydrox- 
ide was obtained from Qualigen Fine Chemicals, India. Freshly 
prepared distilled water was used in all the experiments. 


2.2. Electrode preparation 


The anode electrode was prepared by dispersing the required 
quantity of catalyst powder in Nafion® solution for 30 min using 
an ultrasonic water bath to the prepare catalyst slurry. The cata- 
lyst slurry was spread on a carbon paper in the form of continuous 
wet film. It is then dried in an oven for 30min at 80°C. The 
catalyzed carbon paper was stuck on to the nickel mesh using 
Teflon® dispersion. The prepared electrode was pressed at a 
pressure of 20 kecm "7 and temperature of 120°C for 5 min 
to obtain the composite structure of the electrode. The area of 
the working electrode was 10cm”. Finally, the composite was 
heated at 573 K for 4h to obtain final form of the anode elec- 
trode. The cathode electrode was prepared in a similar manner 
using Monti as catalyst to obtain the composite structure. An 
additional carbon paper was stuck on to air-side of the electrode 
with the help of Teflon dispersion. Finally, a Teflon film was 
stuck on to the carbon paper to obtain the final form of the cath- 
ode electrode. Thus the cathode consist of four layers, i.e., the 
catalyst layer on the carbon paper, nickel mesh, the second layer 
of carbon paper and the Teflon layer. The porous Teflon layer 
on air side of the cathode helps in restricting KOH electrolyte 
transport to the air side and at the same time allows diffusion of 
oxygen from air to catalyst layer. 


2.3. Experimental set-up 


The schematic diagram of the experimental set-up is shown 
in Fig. 1. The prepared anode was placed at the bottom of a 
glass beaker. A hollow cylinder with cathode attached to one 
of the open ends, was placed inside the beaker. Two wires con- 
nected to anode and cathode were acted as terminals. The space 
between anode and cathode was filled with mixture of elec- 
trolyte and fuel solution. The fuel cell was kept in a water bath 
(not shown in Fig. 1) to keep the temperature of the fuel cell 
constant. The fuel-electrolyte and exhausted-fuel storage tanks 
were connected to inlet and outlet of the fuel cell by flexible tubes 
through peristaltic pumps. A magnetic stirrer was used to main- 
tain uniformity of fuel-electrolyte mixture. To protect the anode 
from the stirrer, the anode was shielded with the help of a perfo- 
rated plastic support. The terminals of the cell were connected to 
the digital multimeter (Sanwa PC5000, Japan) through variable 
load. The oxygen present in air acted as oxidant. The detail of 
the experimental set-up and procedure are given in Verma et al. 
[11,12]: 


2.4. Method 


Half-cell analyses of Pt-black anode and MnO?» cathode were 
conducted using cyclic voltammetry. The optimize conditions 


Fig. 1. Schematic diagram of the direct alkaline fuel cell. (1) Fuel-electrolyte 
storage; (2) exhausted-fuel-electrolyte storage; (3 and 4) peristaltic pumps; (5) 
load; (6) anode terminal; (7) cathode terminal; (8) air; (9) cathode electrode; 
(10) anode electrode; (11) fuel and electrolyte mixture; (12) magnetic stirrer; 
and (13) anode shield. 


of fuel and electrolyte concentrations were determined for max- 
imum current density. The alkaline fuel cell constructed was 
fed with mixture of given concentrations of fuel (methanol or 
ethanol or sodium borohydride) and electrolyte (KOH) solution 
from the storage tank at the rate of 1 mL min~!. The fuel and 
electrolyte solution was maintained up to a certain level such 
that the one side of the cathode was in contact with electrolyte 
while the other side exposed to air. The solution was continu- 
ously stirred with a magnetic stirrer, which was located on the 
perforated anode shield. The heat of the reaction is carried out 
of the fuel cell by the un-reacted fuel, products and electrolyte 
mixture. The open circuit voltage (OCV) and the short circuit 
current were measured with the help of a multimeter. The volt- 
age and corresponding current was measured at different loads 
to obtain Lu (current density—cell voltage) characteristic curve. 
The steady-state values of voltage and current were noted in 
all the experiments and the reproducibility of data is checked. 
The influence of different experimental conditions, e.g., elec- 
trolyte concentration, fuel concentration and temperature on fuel 
cell performance were noted. The area resistance of the fuel 
cell was measured by the current interruption technique using a 
potentiostat. 


3. Model development 


Electrical energy is obtained from a fuel cell when a current 
is drawn, but at the same time, cell voltage drops due to vari- 
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ous irreversible losses. The performance of fuel cell is generally 
reduced due to the loss of some electrical energy to drive the oxi- 
dation and reduction reaction at the anode and cathode, which is 
known as activation overpotential, and ohmic resistance of the 
electrolyte, electrolyte—electrode interface and the electrodes. 
Further, the loss occurs due to the mass transfer resistance expe- 
rienced by fuel and oxidant to reach the anode and cathode, 
which is known as concentration overpotential. These losses are 
denoted by nac as the activation overpotential, noh as the ohmic 
overpotential and conc as the concentration overpotential. Thus, 
the cell voltage, Ecen, is written as, 


Ecen = E — (Nac + Noh + Neonc) (1) 


where E is the open circuit voltage. Before the expressions for 
each overpotential losses are derived, the mechanisms of electro 
oxidation and reduction reaction are discussed. 


3.1. Mechanism of electro-oxidation and electro-reduction 


The electro-oxidation reaction mechanism of methanol as 
proposed by Tripković et al. [16] is used to model the activa- 
tion overpotential. Tripković et al. [16] suggested a dual-path 
reaction mechanism (Eq. (2) during electro-oxidation of alco- 
hols (methanol, ethanol, n-propanol and n-butanol) in an alkaline 
solution. 


Alcohol ——» reaction intermediate 


poisoning species (CO) —> CO, 

The end product of methanol electro-oxidation reaction is 
formic acid. However, in alkaline medium formic acid exists in 
the form of formate ion. Apart from formic acid a small quantity 
of carbon dioxide is also formed [16] and it is neglected in the 
model development. In the literature, the report on identification 
of the products formed during oxidation of methanol in alkaline 
solution is scanty. The path shown by Eq. (2a) is the main reac- 
tion path and assumes the formation of reactive intermediates, 
weakly bonded to the surface [17]. 


CHOH & HCO, + 3H* + 3e 


+30H™ 
3 H2O (3) 
OH” = Oldie (4) 
HCOad + OHag > HCOOH (5) 


Eq. (3) is a dehydrogenation process involving several steps 
and may be written as follows [18,19]: 


—H(H*++e7) 


HH +e- 
CH3;0H > HE ten 


—H(Ht+ = 
Hte CH Oa Ee 


—> CHA 
(6) 


Only a few literatures are available on the ethanol electro- 
oxidation. The C-C bond breaking for the complete oxidation of 
ethanol impedes the process of ethanol electro-oxidation. There 
is no electrocatalyst, which enables complete ethanol oxidation 
in the temperatures range of 30-90 °C [15,20,21]. Acetic acid 


CH3 02a 


—-» Acid (in anionic form) 


forms as a final product of ethanol electro-oxidation. A small 
quantity of carbon dioxide generated as reported by Tripković et 
al. [16] and Tarasevich et al. [21] is neglected in the model devel- 
opment. The reaction mechanism for ethanol electro-oxidation 
is similar to that of methanol electro-oxidation [16], which is 
given by, 


C)H,OH © CH + 3H* + Ze 


+30H™ 
Lee. 3 H20 


(7) 
OH" & OHaa+e7 (8) 
CH3COaq + OHag > CH3COOH (9) 


Morris et al. [22] reviewed the electro-oxidation of boron 
compounds and mentioned that the anodic oxidation of BA" 
in alkaline medium on platinum electrode occurs via ioniza- 
tion of hydrogen from the parent ion (BH47 ) and from the first 
hydrolysis product, BH30H , bound to the electrode surface. 
They pointed out similar mechanism of borohydride oxidation 
as that of methanol. Lee et al. [9] studied the characterization 
of an alkaline fuel cell that uses NaBH; as fuel. The reaction 
mechanism suggested by them is given below: 


BH4 os HBOna4 +5H20 + Be (10) 
(2a) 
(2b) 
OH” & OHag +e7 (11) 
HBO2a4 + OHag +e —> BO 2+H20 (12) 


The reduction of oxygen in alkaline medium generally pro- 
ceeds by either of the two pathways [23,24]. They are described 
below as direct oxygen reduction to OH ions, i.e., a four- 
electron pathway: 


O2 +2H20 + 4e° > 40H (13) 


or, an oxygen reduction to HO2~ ions, i.e., two electron 
pathway: 


O2 +H20 + 2e° — HO »+0OH (14) 


with subsequent reduction of peroxide ion to OH” ions, i.e., 
two electron pathway: 


HO 2 +H20 + 2e — 30H (15) 


Eqs. (14) and (15) are collectively produce 2+2 elec- 
tron mechanism [23]. The cyclic voltammmetry and oxygen 
consumption experiments have shown that the oxygen electro- 
reduction reaction follows four-electron transfer mechanism 
[11]. 


3.2. Activation overpotential 


Activation overpotential is directly related to the electro- 
chemical reaction kinetics since the reaction propagates at the 
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rate demanded by the specific load. At no load condition, 
oxidation-reduction reactions reach equilibrium at the elec- 
trodes. The current density related to electron transfer in such 
oxidation—reduction reactions is known as exchange current. 
Many investigators [3—6,25] studied the cathode reaction in alka- 
line condition and reported that the oxygen reduction kinetics is 
more favorable at cathode as compared to the oxidation of fuel 
at anode. Thus, it is assumed that the activation overpotential 
at cathode is less significant compared to that at anode. In the 
following sections, the expressions for activation overpotentials 
are derived for different fuels. 


3.2.1. Methanol 

The electro-oxidation reaction mechanism of methanol (Eqs. 
(3)-(5)) as proposed by Tripković et al. [16] is used to model 
the activation overpotential. Eqs. (3) and (4) are assumed much 
faster than the overall reaction rate-determining step (Eq. (5)). 
In our model development the reaction order with respect to 
methanol concentration is 1.0, whereas reaction order with 
respect to OH is taken as 0.5 [18]. The rate expressions for 
Eqs. (3)-(5) are given below: 


rı = kiCm — by fuco (16) 
r2 = bag — kh Oon (17) 
r3 = k3Oucoin (18) 


In Eqs. (16)-(18), kı, k2 and k3 are the rate constants for 
the forward reactions (Eqs. (3-5)) and E and k are the rate 
constants for the backward reactions (Eqs. (3) and (4)). The 
expression for k3 is given by [26]: 


an Fnac 
RT 


k3 = kK exp ( (19) 
where o is the transfer coefficient, which is defined as the pro- 
portion of the electrical energy used for changing the rate of 
electrochemical reaction (Eq. (5)). Its value depends on the 
nature of the reaction and the electrode material and its value 
must be between 0 and 1.0. The value for œ in the case of 
methanol oxidation is found as 0.22. Wang et al. [25] used a 
value of 0.24 for a similar reaction. n, the number of electrons 
transferred during oxidation of methanol, is 4. Here, K9 is the 
standard rate constant, F the Faraday constant, T the tempera- 
ture and R is the universal gas constant. The concentration of 
methanol and OH™ are represented by Cy and Cop (mol LT!), 
respectively. r; (i= 1-3) is the rate of reaction and 6; (i= OH, 
HCO) is the fractional coverage defined as the ratio of number 
of adsorbed species to the total number of surface sites of the 
catalyst. At equilibrium, rı and r2 is equated to zero and by 
solving the Eqs. (16)—(18) it follows: 


r= k3Kı KoCuCony (20) 
where Kı GE) and K2 (=k2/k5) are adsorption—desorption 


equilibrium constants correspond to Eqs. (3) and (4), respec- 
tively.Also, the net rate of the electrode reaction is governed by 


the rate-determining step hence it is expressed as: 
E 21 
bebe (21) 


where j is the current density. On solving Eqs. (19)-(21), the 
relation between activation overpotential and current density for 
methanol is given by, 


RT jen Cor 
es l 22 
Nac (=) n ( Ku , (22) 


where KM (=nFk?K 1K2) is a constant for methanol oxidation. 


The appropriate values for k, Kı and K3 are required to estimate 
Ky. 


3.2.2. Ethanol 

The reaction mechanism (Eqs. (7)-(9)) proposed by Trip- 
ković et al. [16] is used to model the activation overpotential 
for ethanol electro-oxidation. The reaction order with respect 
to ethanol concentration is taken as 1.0, whereas reaction order 
with respect to OH is taken as 0.5 [21]. The corresponding rate 
expressions for Eqs. (7)—(9) are given below: 


r4 = bar — k4OcH,co (23) 
rs = bet — k; on (24) 
re = k6OcH,C0 0S (25) 


where k4, ks and ke are the rate constants of the forward reactions 
(Eqs. (7)-(9)), k4 and ks are the rate constants of the backward 
reactions (Eqs. (7) and (8)). Following similar method that car- 
ried out in the case of methanol and using o value of 0.3 and 
n=4, the relation between activation overpotential and current 
density for ethanol is given by, 


RT Je Con 
a ji : 26 
Nac E =) n ( Ke (26) 


where Kr Gab: k? K4K5) is a constant for ethanol oxidation. 


3.2.3. Sodium borohydride 

The reaction mechanism (Eqs. (10)-(12)) suggested by Mor- 
ris et al. [22] is taken into consideration and the corresponding 
rate expressions are given below: 


rı = kıCspg — ky Oupo, (27) 
rg = bet — kk Oon (28) 
ro = ofge, Cé (29) 


where k7, kg and ko are the rate constants of the forward reac- 
tions (Eqs. (10)-(12)), k4 and kg are the rate constants of the 
backward reactions (Eqs. (10) and (11)). The order of the reac- 
tion with respect to borohydride is taken as 1.0 and 0.5 with 
respect to OH ion [22]. Using similar procedure as described 
before and taking «œ =0.3 and n = 8, the relation between activa- 
tion overpotential and current density for sodium borohydride 
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is given by, 


RT l 
= SCC n 
Nac DES 


where Ksg (=nF. HK: Kai is a constant for sodium borohydride. 


yaad. A055 
JCsp Con 


; (30) 
Ksp 


3.3. Ohmic overpotential 


In modeling of ohmic overpotential, it is assumed that the 
ohmic drop in current collectors and electric connections is 
negligible and change in resistance of electrolyte due to the 
presence of fuel is insignificant. The first assumption is rea- 
sonable since modeling of single fuel cell is involved and stack 
is not considered. The second assumption implies that the ohmic 
losses are same for all the three fuels used. The ohmic losses are 
proportional to the current density and represented as: 


Noh = JRoh, (31) 


where Ron is the area specific resistance offered by electrolyte 
for ion transport. The variation of resistance with the KOH con- 
centration is reported by Jo and Yi [27]. The expression for area 
specific resistance, BA, at 25°C is determined from the regres- 
sion analyses of the experimental data on resistance and KOH 
concentration presented by Jo and Yi [27]. The expression for 
area specific resistance, Rio at 25 °C is found out as follows: 


Roy = 13.1737 — 90.336con + 230.5735chy — 81.07083ch4 
(32) 


It is well known that the resistance to ionic mobility decreases 
in the electrolyte solution with the increase in temperature. The 
final expression modified for area specific resistance in terms 
of KOH concentration and temperature is determined from the 
experimental data and is given by: 


Ron = Ri, — 25.78 + 1.18757; — 0.0062577, (33) 


where coy is the concentration of KOH in g cc7™!, Tı the tem- 
perature in °C and RA is the area specific resistance? in Q em. 


3.4. Concentration overpotential 


Normally, oxygen at the cathode of a fuel cell is supplied in 
the form of air. The oxygen consumption at cathode leads to 
the decrease in its partial pressure at the cathode. The extent of 
decrease in partial pressure will depend on the current drawn 
from the fuel cell. The decrease in partial pressure or con- 
centration of oxygen results in concentration overpotential as 
the oxygen will not be replenished immediately at the cathode 
because of the mass transfer resistance. The similar phenom- 
ena occur in the anode compartment where the oxidation of fuel 


2 Note that unit of length is in cm because the original experimental data 
published by Jo and Yi [25] is in cm and the regression analyses expression is 
valid for the same unit. In our final calculation of ohmic loss, the length unit is 
converted to m. 


Teflon film 


GDL 
Ni-mesh 
GDL 

MnO, layer 


Electrolyte (KOH) + Fuel 


Fig. 2. Schematic diagram of cathode components. 


takes place and fuel concentration decreases. In the experimen- 
tal set-up it is shown that the fuel and electrolyte mixture is 
well stirred. Thus, the effect of concentration overpotential at 
anode is reduced considerably and it is assumed to be insignif- 
icant compared to that at cathode. The schematic diagram of 
the cathode is shown in Fig. 2. It is seen in Fig. 2 that oxygen 
present in air has to diffuse through porous Teflon layer, gas dif- 
fusion layer, Ni-mesh and another gas diffusion layer to reach 
cathode catalyst (MnOz) surface. In estimating the mass transfer 
resistance the thickness of the composite cathode made of gas 
diffusion layers and Teflon layer is considered. It is assumed that 
the Ni-mesh offers negligible resistance. It is further assumed 
that the concentration of oxygen at the interface of air and gas 
diffusion layer is Co, i.e. the bulk concentration of oxygen in the 
air, while the concentration of oxygen at the catalyst layer is C. 
Under steady-state condition and from species conservation, the 
oxygen molar flux, No, is related to the current density, j, by: 
J 

No = 75 (34) 

One-dimensional Fick’s law in terms of the molar flux of 
oxygen is used to describe the transport of oxygen in the gas 
diffusion layer. It is given by, 

Co — C 
No = — Deft —— 
ô 

where ô is the thickness of electrode, Deg is the effective diffu- 
sivity of the oxygen-nitrogen gas pair in the porous media with 
the porosity £. Deft is evaluated from the bulk diffusivity, Do air, 


using Bruggeman’s correction [28]: 


(35) 


Dest = €7/* Do air (36) 


The catalyst layer is regarded as a thin film located on the gas 
diffusion layer. The rate of the electrochemical reaction within 
the catalyst layer is described using Butler- Volmer equation [26]. 
It is then simplified in terms of oxygen concentration at the 
cathode catalyst layer as [24]: 


win C exp Dn oun F 
J= Jo Co RT 
where £ is the transfer coefficient and jo is the exchange cur- 
rent density. The value of ĝ varies from O to 1.0 and in the 
present case it is assumed as 0.3. The value of n depends on 


the stoicheometry of oxygen reduction reaction. Since cathode 
reaction (O2 +2H,0+4e  =40H ) is same for all the fuels, 


(37) 
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i.e., methanol, ethanol and sodium borohydride, the value of 
is 4. The exchange current density is given by [25]: 


73200 / 1 1 
) = 0.0422 , 38 
Jo exp | R (= sl (38) 


where Tı is temperature in °C. On solving Eqs. (34)-(37), the 
expression for concentration overpotential is given by, 


RT J e 
anF | pO- jM) GY 


where M is a constant as defined in Eq. (40) and A is the area of 
the electrode. 
AôCo 


M = —— 40 
n Fe?/2 Do air o 


Neonc = 


3.5. Model equations 


In the model development, it is assumed that the direct alka- 
line fuel cell is operated at steady state and isothermal condition 
is prevailed. The final reaction for the direct methanol alkaline 
fuel cell is given by the adding two half-cell reactions taking 
place at anode (Eq. (3)-(5)) and cathode (Eq. (13)) and it is 
given by, 


CH30H + O2-— HCOOH + H20 (41) 


By incorporating the expressions for activation (Eq. (22)), 
ohmic (Eq. (33)) and concentration (Eq. (37)) overpotentials in 
equation (Eq. (1)), the final model equation is obtained for direct 
methanol alkaline fuel cell as, 


RT JOVE CN , 
Ecen = E ( ) ln ( M “OH Jan 


on bh K 


RT ( j ) 
In | - (42) 
anF Jo — jM) 


Similarly by adding the half-cell reactions of ethanol (Eq. 
(7)-(9) and sodium borohydride (Eq. (10)-(12)) to oxygen 
reduction reaction (Eq. (13)), the reactions for direct ethanol 
and borohydride alkaline fuel cells are given by, 


C.Hs50H + O2 —> CH3COOH + H20 (43) 
BH4 +202— BO? +2H20 (44) 


The model equations for ethanol and sodium borohydride are 
obtained from Eqs. (26), (33), (37) and (30), (33), (37), respec- 
tively. The model equation for direct ethanol alkaline fuel cell 
is given by: 


RT iC. Ca 
Ecen = E ( ) In ( E OH JRoh 


on bh KI 


Bl) ( j ) (45) 
anF "Lat jM) 


The model equation for direct borohydride alkaline fuel cell 
is given by: 


RT jc. Ca ; 
En = E — (=) In Le — JRon 


ued l j (46) 
n 
anF jol — jM 


From the above equations (Eqs. (42), (45) and (46)), the 
operating voltage and current density at a given load of the 
direct alcohol or borohydride alkaline fuel cell is determined 
at the known rate constants and operating conditions, e.g., fuel 
concentration, electrolyte concentration and temperature. 


4. Results and discussion 


The experimental data on cell voltage and current density at 
different operating conditions are presented and compared with 
model predictions. 


4.1. Electrolyte concentration 


Fig. 3a—c shows the experimental data and model predictions 
on current density—cell voltage (j—v) characteristics at four dif- 
ferent KOH concentrations for methanol, ethanol and sodium 
borohydride at 25 °C. It is seen in the figures that the cell volt- 
age increases with the increase in KOH concentration in the 
range of 1-3 M and then it decreases with the further increase 
in KOH concentration. The reason for the decrease in voltage 
may be the less availability of methanol at the catalyst site with 
the increase in KOH concentration. The experimental data on 
Lu plot shown by symbols are predicted by the model as shown 
by the lines for different KOH concentrations. The model pre- 
diction in the region of concentration overpotential shows larger 
deviation compared to the regions of activation and ohmic over- 
potentials. It may be because of the assumption made in the 
model that the concentration overpotential at anode is negligi- 
ble compared to that at cathode. The abrupt dip in the cell voltage 
at a higher current density is not predicted by the model because 
of the above assumption and perhaps due to undesired methanol 
oxidation at cathode at a high current density. 


4.2. Fuel concentration 


Fig. 4a—c shows current density—voltage characteristics, 
where the experimental data are shown by symbols and the 
model predictions are shown by lines, for four different fuel 
concentrations of methanol, ethanol and sodium borohydride, 
respectively. Itis seen in the figures that the cell voltage increases 
with the increase in fuel concentration for a given value of cur- 
rent density. Although the cell performance increases initially 
but it does not increase proportionally with further increase in 
fuel concentration. The reasons may be because of the decrease 
in hydroxyl ion mobility with the increase in fuel concentra- 
tion. Also, the active sites for hydroxyl ion adsorption might 
have been blocked by the fuel. Although 3 M alcohol concentra- 
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Fig. 3. Comparison of experimental data with model predictions pertaining to 
current-voltage characteristics for 2 M fuel at different electrolyte (KOH) con- 
centrations (T; =25 °C): (a) methanol (b) ethanol and (c) sodium borohydride. 


tion generated slightly higher cell performance, the difference in 
performance between 3 and 2 M alcohol concentration is negligi- 
ble. Subsequently, experiments were carried out at 2 M alcohol 
concentration to study the effect of temperature. The current 
density-cell voltage plots are somewhat predicted by the model 
for different fuel concentrations. The model overpredicts the 
experimental data in the region of higher current density, where 
concentration overpotential is dominant. The higher current den- 
sity originates because of the lower load and consequently faster 
reaction kinetics expected at the anode and cathode. The sharp 
decrease in cell voltage value at a higher current density com- 
pared to that predicted by the model may be because of the oxida- 
tion of fuel at cathode at a higher fuel concentration, especially 
when system demands higher reaction kinetics at the electrodes. 


4.3. Temperature 


Fig. 5a—c shows the experimental data and model predictions 
on current density—voltage characteristics at three different oper- 


ating temperatures for 3M KOH and 2M methanol, ethanol, 
sodium borohydride, respectively. It is seen in the figures that 
the cell voltage increases with the increase in temperature for a 
given value of current density. The performance increases with 
the increase in temperature because of decrease in activation 
overpotential, concentration overpotential and increase in ionic 
conductivity and mobility at a higher temperature. The perfor- 
mance of direct sodium borohydride alkaline fuel cell does not 
increase appreciably with the increase in temperature and in 
fact shows decreasing trend in current—voltage characteristics at 
65 °C. The reason for this decrease may be because of hydrogen 
gas generation from sodium borohydride at higher temperature. 
The hydrogen gas generated escapes from the cell without react- 
ing. The maximum power density of 24.3 and 14.5 mW cm~? 
were obtained for methanol and ethanol at 65°C while that for 
NaBH; was 22.5 mW cm~? at 25°C. The experimental data, 
shown by the symbols, are reasonably predicted by the model as 
shown by the lines except at a higher current density. The reason 
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Fig. 4. Comparison of experimental data with model predictions pertaining 
to current-voltage characteristics at different fuel concentrations (T; =25 °C; 
Con =3 M): (a) methanol (b) ethanol and (c) sodium borohydride. 
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Fig. 5. Comparison of experimental data with model predictions pertain- 
ing to current-voltage characteristics for 2M fuel at different temperatures 
(Cou =3 M): (a) methanol (b) ethanol and (c) sodium borohydride. 


for overprediction of experimental data at a higher temperature 
is already discussed. 

A consistent overprediction of experimental data by the 
model at a higher current density has led to the idea of addi- 
tion of an excess concentration overpotential term to take into 
account of various undesired losses at the electrodes, e.g., oxida- 
tion of fuel at cathode and concentration overpotential at anode. 
The additional term used in the model is given in the appendix. 
Fig. 6a—c shows that the modified model predicts the experi- 
mental data well in the range of higher current density values 
for methanol and ethanol. Fig. 6a and b shows the prediction by 
the modified model for methanol at different KOH and fuel con- 
centrations, respectively. Whereas, Fig. 6c shows the prediction 
by the modified model at different operating temperatures for 
direct ethanol alkaline fuel cell. It should be noted that source 
of error in modeling may be from the use of Eq. (1) in which 


overpotential terms are added in linear manner. Since kinetics of 
the reaction at anode and cathode depends on concentration, an 
integrated approach should be taken to account for all the losses. 
A consistent over prediction of the cell voltage at a higher cur- 
rent density may have resulted due to assumed linear behavior 
of the loss terms. 


4.4. Model parameters 


The parameters used in the model Eqs. (42), (45) and (46) 
for different fuels, e.g., methanol, ethanol and sodium borohy- 
dride are presented in Table 1. It is seen that the fitted value 
for the adsorption and desorption constants, Kı and K3 of Eq. 
(20) and the rate constant, H (Eq. (19)), of the rate determining 
step is of the same order of magnitude to that found in literature 
for similar reaction involving methanol at 70°C. The constant, 
Kg, for ethanol as listed in Table 1 is of the same order that for 
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Fig. 6. Comparison of experimental data with predictions by the modified model 
pertaining to current-voltage characteristics (Tı = 25 °C unless otherwise speci- 
fied): (a) 2 M methanol at different KOH concentrations (b) 3 M KOH at different 
methanol concentrations (c) 2 M ethanol and 3 M KOH concentrations at differ- 
ent temperatures. 
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Table 1 
Model Parameters 
Parameter Temperature Unit Present model value/literature value Reference 
Kı 65°C mol`! L 1.0/1.0 (at 70°C) Sundmacher et al. [17] 
Ko 65°C mol IL 0.01/0.01 (at 70°C) Sundmacher et al. [17] 
H 65°C mol m~? s7! 0.155/0.01 (at 70°C) Sundmacher et al. [17] 
Km (=nFK, K2 ko) 25°C mol`! L 3.0/not available Sundmacher et al. [17] 
45°C 5.0/not available 
65°C 6.0/38.6 (at 70°C) 
Kg (=nFK4K5 ke) 25°C mol`“! L 1.0/not available Not available 
45°C 2.3/not available 
65°C 4.0/not available 
Ksg Gank kat) 25°C mol m~? s7! 0.001/not available Not available 
45°C 0.03/not available 
65°C 0.25/not available 
a Dimensionless 0.22/0.24 Wang et al. [23] 
a! Dimensionless 0.3/0.5 Tarasevich et al. [19] 
a” Dimensionless 0.3/not available Not available 
B Dimensionless 0.3/0.1-0.5 Larminie and Dicks [27] 


methanol, Kym. It is expected that the rate constants of electroox- 
idation reaction and adsorption and desorption constants for two 
consecutive alcohols in the homologous series should be of the 
same order of magnitude. Since the rate constants for eletro- 
oxidation and adsorption and desorption constants for ethanol 
and NaBH; on electrode surface are not reported in the litera- 
ture, they are not compared in Table 1. The «œ values for different 
fuels and the $ value used are within the comparable range as 
quoted in the different literatures [19,25,29]. 


5. Conclusions 


The direct alkaline fuel cell constructed for the use of 
methanol, ethanol, sodium borohydride as fuel was tested and 
current-voltage characteristic curves were obtained at differ- 
ent experimental conditions, e.g., electrolyte concentration, fuel 
concentration and temperature. Pt-black was used as the anode 
catalyst and MnO> was used as the cathode catalyst. KOH solu- 
tion was used as the electrolyte. The cell performance increases 
initially with the increase in electrolyte (KOH) concentration and 
then decreases with the further increase of the same. Initially, the 
cell performance increases with the increase in fuel concentra- 
tion, but no appreciable increase in the performance is recorded 
at a higher fuel concentration. The performance of the fuel cell 
increases with increase in temperature in general with the excep- 
tion to NaBH; alkaline fuel cell. A mathematical model for the 
prediction of voltage at a given current density of direct alkaline 
fuel cell is developed, by taking into account the losses due to 
activation, ohmic and concentration overpotentials. The activa- 
tion overpotential term is formulated using reaction mechanisms 
proposed earlier in the half-cell studies. Ohmic overpotential is 
modeled based on conductivity of electrolyte solution (KOH) 
at different temperatures. The concentration overpotential is 
formulated using Fick’s law and Butler-Volmer equation. The 
model reasonably predicts the experimental data on cell voltage 
and current density and captures the influence of electrolyte and 


fuel concentrations and temperature with the exception to data 
at a higher current density. Addition of an excess concentration 
overpotential term to the model equation further improves the 
model prediction for the given range of current density. 
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Appendix A. Undesired losses at the electrodes 


The direct alkaline fuel cell has suffered from minute losses at 
the electrodes in the range of high current density. These losses 
are concentration overpotential at anode and undesired fuel oxi- 
dation at cathode as stated earlier. An important observation in 
the modeling study is that the difference between model pre- 
diction and experimental data increases with the increase in the 
current density. At a higher current density, the fuel at the anode 
is used up at a faster rate not allowing to be replenished because 
of the mass transfer resistance. An excess loss term as a function 
of current density is added to the model equation to quantify the 
various losses [30]. This is given by, 


nFL = L x exp(N x j) (A.1.1) 


where up is the excess concentration overpotential, L 
and N are functions of the electrolyte concentration. L 
and N have been obtained empirically for the different 
fuels from experimental data. The constants L and N for 
methanol are found out by regression analyses as, L = 1 x 
1077 x e!45Con, N = 0.0026C2,, — 0.025Coxn + 0.0877 and 
that for ethanol, L = 9 x 10-7 x e23CoH, N = 0.0021C2), — 
0.025CoH + 0.0758 and for sodium borohydride, L = 2 x 
107° x el 5Con, N = 0.0014Cĝy — 0.015Con + 0.079. The 
regression coefficient is 0.985. 
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